1, 2008 20:17 WSPC/Trim Size: 9in x 6in for Proceedings 



calll 



ANALYTIC CALCULATION OF 
BOSE-EINSTEIN CORRELATIONS AT LEP1 AND LEP2 * 



T. CSORGO 

MTA KFKI RMKI, H-1525 Budapest 114, P.O.Box 49, Hungary 
E-mail: csorgo@sunserv.kfki.hu 

Theoretical results in the r-model are discussed in the context of Bose-Einstein 
correlations in e + e — annihilations at LEP1 and LEP2 as well as that of relativistic 
heavy ion collisions at SPS and RHIC. 

A current motivation 1 to study of Bose-Einstein 2 (or HBT 3 or GGLP 4 ) 
correlations in e + e~ annihilations at LEP1 or LEP2 is to determine their 
possible effects on the mass reconstruction of W bosons at LEP2, see e.g. 
ref. 5 for a recent review. The goal is to progress in 2 steps 5 : 

1) understand Bose-Einstein correlations at LEP1, at the Z° peak, 

2) generalize these results to the case of W + W~ production at LEP2. 

1. Bose-Einstein correlations in "r-model at LEP1 

A model of strongly correlated phase-space was developed in ref 6 to explain 
the experimentally found invariant relative momentum Q = y/—(ki — fe) 2 
dependence of Bose-Einstein correlations (BEC-s) in e + e~ reactions 7 ' 8 . As 
the correlation function was expressed in terms of a distribution of a proper- 
time t, let us call this model the r-model 6,9 ' 10 . It is based on 

Assumption i): The momentum k and the average position x are 
strongly correlated : x = dk where d is a constant of proportionality, 
specified below, and x = x^ L — (t, r), k = k^ — (E^, k). 

Assumption it): The space-time momentum-space correlation is nar- 
rower than the proper-time distribution (when both are measured in di- 
mensionless units). 
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science fellowship and by the MTA-OTKA-NSF grant INT0098462. 
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Thus the emission function of the r-model is 

/>oo 

S(x,k) = H(t)6 a (x- dfc)JVi(fc) 
Jo 

where H(t) is the proper-time distribution, the factor 5 a {x — dk) is de- 
scribing the strength of the correlations between coordinate space and mo- 
mentum space variables. In the original version of the T-modcl, these cor- 
relations were taken to be ideally strong, 5a (x — dk) = 5^(x — dk), a 
four-dimensional Dirac-delta. The results of the r-model are unchanged 9 
if the (x, k) correlations are described by a sufficiently, but not infinitely 
narrow function, whose width is characterized by some dimensionless and 
narrow scale A <C At/t, where At is the width of H(t). The distributions 
are normalized, J °° H(t) = 1, and J d 4 x5A(x — rk/m) = 1. 

If the expansion is 1+1 dimensional, then d = Ti/mt, if the expansion 
is 1+3 dimensional then d = rjm, where the (longitudinal) proper-time is 
defined as (t; = yjt 2 — r 2 and) r = ^Jt 2 — r 2 — r 2 — rf , while the (trans- 
verse) mass is given as (m t = \J E\ — k. 2 and) m = y/E^. — k 2 . The former 
is relevant in 2-jet decays of Z° at LEP1, the latter is relevant in the case 
of fully hadronic W + W~ decays, mostly 4-jet events. 

The experimentally measurable single-particle spectra N\{k) is 

J d i xS(x,k) = JVi(fc) 

arbitrary, and can be taken directly from the measurements. 

The two-particle BEC-s C 2 {k 1 ,k 2 ) = N 2 {ki,k 2 )/(N 1 (ki)N 1 (k 2 )) are 
calculated with the help of the Yano-Koonin formula 6 ' 9 ' 10 . The key step is 
that for any choice of d, one gets (xi — x 2 )(k\ — k 2 ) = — 0.5(d\ + d 2 )Q 2 . 
Using a saddle-point integration based on i) and it), we find that 

C 2 (h,k 2 ) ~ 1 + XRe[H 2 (w)], 

where the argument is w — Q 2 /(2m) for 1+3 dimensional, and w = 
Q 2 /(2m t ) for 1+1 dimensional expansions, corresponding to the appro- 
priate choice for variable d, and H(w) — J 00 dTH(r) exp(iwr) is a Fourier- 
transform of H(t). Thus an invariant relative momentum Q dependent 
BEC appears, if the (x, k) correlations between spacetime and momentum- 
space are strong enough and the proper-time distribution is broad enough. 

For example, let us consider two-jet events at LEP1. These are 1+1 
dimensionally expanding systems, hence d = Ti/m t - In case of the asym- 
metric Levy H(t) distribution, H(t) oc ,_\ 3 /2 exp(— ), we get 6 

C(Q) = 1 + cos(R 2 Q 2 ) exp{-QR) 
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with Rq = To/mt and R 2 = Ar/mt- This form is similar to the observed 
data at LEP1. Currently known 5 experimental constraints on BEC-s in 
e + e" annihilations at LEP1 are summarized below. It is easy to show that 
the r-model satisfies all these properties 9 ' 10 : 

(1) BEC-s exist in both Z° decays at LEP1 and in fully hadronic decays 
of W+W- events at LEP2. 

(2) The pion source is not spherically symmetric, but elongated along 
the thrust axis, when analyzed in the LCMS frame. 

(3) The shape of the BEC function is far from Gaussian at LEP1. 

(4) The effective source size decreases as 1/ ^fm for pions, kaons, protons 
and A particles, and as 1/^/rnt for pions with different m*. 

(5) The normalized three-particle cumulant correlation coincides with 
1 within experimental errors in a large region of Q3. 

(6) No significant isospin dependence is seen. 

(7) In e+ + e~ annihilations, the two-particle BEC depends on the rela- 
tive momentum components only through the invariant momentum 
difference Q — ^—{k\ — k 2 ) 2 - 

(8) After the peak, the BEC develops a shallow C(Q) < 1 dip region. 

2. Bose-Einstein correlations from the r-model at LEP2 

The key step is to introduce r for e+e~ — ► W + W~ fully hadronic decays. 
This variable is now measured from the production point of the W + W~ 
pair (and not from the decay point of any of the W-bosons). 

In these reactions the particle emitting source has two components, 

S(x,k) = S + (x,k) + S-(x,k), 

corresponding to the decays of W + and W~ , respectively. Thus the sin- 
gle particle spectrum has also two components, N\(k) — N^(k) + N^(k). 
When calculating the BEC-s for such a binary source, the key observation is 
H + (t) = H~{t) = H 2 (t). This simple equation has dramatic consequences 
for the BEC of fully hadronic decays of W + W~ pairs at LEP2: 

C 2 {k u k 2 ) = l + \Re[H 2 2 (^-% 

where H 2 {w) = / °° dTH 2 {T) exp(iwr) = if ± (w) is the Fourier-transformed 
proper-time distribution at LEP2. If H 2 {t) ~ H\{t), where the latter 
stands for the proper-time distribution of particle production in e + e~ — > Z° 
hadronic decays at LEP1, then the BEC-s at LEP1 and LEP2 are also 
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approximately the same. This result implies no genuine inter-W Bose- 
Einstein correlations within the r-model, but the existence of a fully chaotic 
source with full symmetrization, see ref 10 for further discussions. 

3. Implications for heavy ion collisions 

The BEC-s or HBT correlations at SPS and at RHIC were mostly analysed 
in terms of the Bertsch-Pratt (BP) or side-out-long variables, see refs 11 ' 12 
for recent reviews. If we assume that the r-model and strong (x, k) corre- 
lations are relevant also in heavy ion collisions, we can decompose the re- 
sulting correlation functions in the BP variables to find a surprizing result: 
The variable of the correlation function, Q 2 Ar/m t = Q 2 R 2 + Q 2 R 2 + Q 2 R 2 
yields R 2 = Ar/m t , Rf = Ar/mt and R 2 = m 2 Ar/mj, where subscripts 
{ Sl0 ,i} stand for side, out and long, respectively. 

Experimentally, this implies R s ss Ri oc \/^Jm~t, and R decreases with 
m t faster than these. Thus very strong (x, k) correlations in the r-model 
generate R < R s dynamically! More detailed studies are needed to deter- 
mine the proper variables of BEC-s in high energy heavy ion collisions. If 
the r-model is relevant in high energy heavy ion collisions, it may imply 
that R <C R s means a very broad proper-time distribution and strong (x, k) 
correlations. Experimentally, this question can be decided by determining 
if the Bose-Einstein correlation function at RHIC depends on the relative 
momentum components only through the invariant relative momentum Q, 
as in the case of e+e~ collisions, or not. 
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